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The study of d(d,p)t reaction is very important for the nucleosynthesis in both standard Big Bang 
and stellar evolution, as well as for the future fusion reactors planning of energy production. The 
d(d,p)t bare nucleus astrophysical S(E) factor has been measured indirectly at energies from about 
400 keV down to several keV by means of the Trojan horse method applied to the quasi-free process 
2H (SLi, pt)*He induced at a lithium beam energy of 9.5 MeV, which is closer to the zero quasi-free 
energy point. An accurate analysis leads to the determination of the Spare(0) = 56.7 +2.0keV -b and 
of the corresponding electron screening potential Ue = 13.2+4.3eV. In addition, this work gives an 
updated test for the Trojan horse nucleus invariance comparing with previous indirect investigations 
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using ?He = (d + p) breakup. 


PACS numbers: 26.20.Cd, 25.45.Hi, 24.50.--g 


I. INTRODUCTION 


The d+d nuclear reactions are important in both nu- 
clear astrophysics and fusion energy applications 

These reactions are among the thermonuclear pro- 
cesses occurring during the first minutes of the universe 
immediately after the Big Bang. In particular, knowl- 
edge and modelling of the primordial abundance of deu- 
terium, which depends on precise cross section data, give 
important information about the baryon density of the 
universe. Moreover, primordial deuterium is burned dur- 
ing the earliest evolution stage of stars: the pre-main 
sequence phase. Thus, a better knowledge of the pa- 
rameters characterizing these reactions can improve our 
understanding of the first phases of stellar evolution. As 
for the Standard Big Bang Nucleosynthesis, the region of 
interest ranges from 50 to 300 keV, and it is only several 
to 20 keV for stellar evolution processes. 


In addition to these important astrophysical topics, 
the interest of scientists around reactions involving deu- 
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terium has been also triggered by the promising possibil- 
ity of exploit them as a powerful and low-polluting source 
of energy in fusion reactors. These reactions belong to 
the network of processes inside the fusion reactors. These 
reactors are expected to operate in the temperature range 
of kT = 1-30 keV. 


Several experiments have been performed below 200 
keV, but available data are not always in agreement 
within each other and some of them are affected by large 
systematic errors. Another weak point is that available 
data below 10 keV, region of interest for fusion reactors as 
well as for burning deuteron in the Pre-Main Sequence 
phase of stellar evolution, are affected by the electron 
screening. 


For these reasons, new indirect experimental studies 
were called for to provide new data in the full range of 
interest for pure and applied physics. The Trojan Horse 
Method (THM) lo, [7] has been applied to the indirect 
study of the d+d reactions using ?He = (d + p) and 
6Li = (d +a) breakup [19,[16], but the °Li breakup data 
give much less points and larger errors than that in the 
case of ?He breakup. 


In this paper, we report on a new investigation of the 
d(d,p)t reaction by means of the THM applied to the 
?H(°Li, pt)^He quasi free process with a beam energy of 
9.5 MeV, which is closer to the zero quasi-free energy 
point. 
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FIG. 1: (Color online) The schematic representation of 
Trojan-horse method 


II. TROJAN HORSE METHOD 


'The Coulomb barrier and electron screening cause dif- 
ficulties in directly measuring nuclear reaction cross sec- 
tions of charged particles at astrophysical energies. To 
overcome these difficulties, the THM Iri has been intro- 
duced as a powerful indirect tool in experimental nuclear 
astrophysics [8121]. The THM provides a valid alter- 
native approach to measure unscreened low-energy cross 
sections of charged particle reactions. It can also be used 
to retrieve information on the electron screening poten- 
tial when ultra-low energy direct measurements are avail- 
able. 

'The basic assumptions of the THM theory have been 
discussed extensively in [2], and the detailed theoret- 
ical derivation of the formalism employed can be found 
in 

A schematic representation of the process underlying 
the THM is shown in Figure [I] The method is based 
on the quasi-free (QF) reaction mechanism, which allows 
one to derive indirectly the cross section of a two-body 
reaction Eq.(1] 


A+z>C+e (1) 


from the measurement of a suitable three-body process 
under the quasi-free kinematic conditions: 


Ata+Ct+er+b (2) 


where the nucleus a is considered to be dominantly com- 
posed of clusters x and b (a = (x & b)). 

After the breakup of nucleus a due to the interaction 
with nucleus A, the two-body reaction (Eq.(I)) occurs 
only between nucleus A and the transferred particle x 
whereas the other cluster b behaves as a spectator to the 
virtual two-body reaction during the quasi-free process. 
The energy in the entrance channel Ea, is chosen above 
the height of the Coulomb barrier ES;P: , so as to avoid 
the reduction in cross section. 

At the same time, the effective energy EA; of the re- 
action between A and z can be relatively small, mainly 
because the energy Ea, is partially used to overcome the 


binding energy e, of x inside a, even if particle x is almost 
at rest the extra-energy is compensated by the binding 
energy of a (Eq.(3)), and the Fermi motion of z inside a, 
Exp, is used to span the region of interest around EZ. 
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Since the transferred particle x is hidden inside the nu- 
cleus a (so called Trojan-horse nucleus), it can be brought 
into the nuclear interaction region to induce the two-body 
reaction A+, which is free of Coulomb suppression and, 
at the same time, not affected by electron screening ef- 
fects. 

'Thus, the two-body cross section of interesting can be 
extracted from the measured quasi free three-body reac- 
tion inverting the following relation: 


do do TE 


— Kw? 5 
dEccdQ pp dC Nm dQ (5) 


where KF is the kinematical factor, |W]? is the momen- 
tum distribution of the spectator b inside the Trojan- 
horse nuclei a, and da/ dQ™ is the half-off-energy-shell 
(HOES) cross section of the two-body reaction A + x > 
C+C: 


do TH do 
A 0C Bug (At > Co) (6) 
l 
where ga (Ax — Cc) is the real on-energy-shell cross sec- 


tion of the two-body reaction A +x — C+ c for the l 
partial wave, P; is the penetration function caused by the 
Coulomb wave function, and Cj is the scaling factor. 


III. EXPERIMENT 


The measurement of the ?H(9Li, pt)*He reaction was 
performed at the Beijing National Tandem Accelerator 
Laboratory at China Institute of Atomic Energy. The 
experimental setup was installed in the nuclear reaction 
chamber at the R60 beam line terminal as shown in Fig- 
ure] The °Li?* beam at 9.5 MeV provided by the HI-13 
tandem accelerator was used to bombard a deuterated 
polyethylene target CDs. The thickness of the target is 
about 160ug/cm?. In order to reduce the angle uncer- 
tainty coming from the large beam spot, a linear target 
with 1 mm width was used. 

A position sensitive detector PSD; was placed at 
40° + 5° to the beam line direction and about 238 mm 
from the target to detect the outgoing particle triton (t), 
and another detector PSD2 was used at 78° + 5? in the 
other side of the beam line at 245 mm from the target 
to detect the outgoing particle proton (p). The arrange- 
ment of the experimental setup was modelled in Monte 
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FIG. 2: Experiment setup of the 7H(°Li, pt)*He reaction 


Carlo simulation in order to cover a region of quasi-free 
angle pairs. A PSD,, was placed at 32° + 5° opposite to 
PSD, as a monitor. The energy resolution of the PSDs 
is about 0.6%-0.8% for 5.48 MeV a source. 

It should be mentioned that no AE detector was 
mounted before PSDs, that was, no particle identification 
was performed in the experiment. It was beneficial to im- 
prove the energy and angular resolution without AE de- 
tectors, because it would lead to additional straggling of 
energies and angles when particles passing through these 
detectors. While it was easy to select events for this reac- 
tion without particle identification from the kinematics 
in data analysis with the help of simulation. 

The trigger for the event acquisition was given by coin- 
cidence of signals by Gate = PSD; x (PSD2+PSD,,). En- 
ergy and position signals for the detected particles were 
processed by standard electronics and sent to the acqui- 
sition system MIDAS for on-line monitoring and data 
storage for off-line analysis. In order to perform position 
calibration, a grid with a number of equally spaced slits 
was placed in front of each PSD for calibration runs. 


IV. DATA ANALYSIS AND RESULTS 


The position and energy calibration of the detectors 
were performed using elastic scatterings on different tar- 
gets (19" Au, !?C, and CD2) induced by a proton beam 
at energies of 6, 7, 8 MeV. A standard a source of 5.48 
MeV was also used. 

After the calibration of the detectors, the energy and 
momentum of the third undetected particle (o) were re- 
constructed from the complete kinematics of the three- 
body reaction Li +d — t+p-+a, under the assumption 
that the first particle is a triton (detected by PSD,) and 
the second one is a proton (detected by PSD»). 


A. Selection of the three body reaction events 


The basic step of data analysis is to select the three- 
body reaction events of 7H(°Li, pt)^He from all exit chan- 
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FIG. 3: (Color online) Selection of the three body reaction 
events of ?H(SLi, pt)*He from the E; — E» kinematic locus 


nels. Figure [3] shows the experimental spectrum of the 
E; — E2 kinematic locus. Comparing with the Monte 
Carlo simulation [1], we can select the range by a graph- 
ical cut (red line polygon) where the three body reaction 
?H(SLi, pt)*He events are located. It will be used as a 
basic selection cut in the following data analysis. 


B. Qs value spectum 


Once selected the three-body reaction events of 
?H(°Li, pt)*He, the experimental Qs value can be ex- 
tracted, as reported in Figure Æ 

There is a peak whose centroid is at about 2.5 MeV 
(in good agreement with the theoretical prediction, Q = 
2.558 MeV). It is a clear signature of the good calibration 
of detectors as well as of the correct identification of the 
reaction channel. 

The events outside of the 2.5 MeV peak came from 
the background and some reaction located in the energy 
range in Figure[3] but not in agreement with the assump- 
tion that the first particle is a triton and the second one 
is a proton, so that the calculated Q3 value deviated from 
the expected value. 

Only events inside the 2.5 MeV Q-value peak were con- 
sidered for the further analysis. 


C. Momentum distribution of a inside °Li 


As in all standard THM analysis, the next step is to 
identify and separate the quasifree mechanism from all 
the other processes. This is usually done by recalling the 
definition of a QF reaction, i.e., a reaction where the third 
particle (spectator) retains the same momentum it had 
within the Trojan horse nucleus. Thus, the momentum 
distribution of the third and undetected particle will be 
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FIG. 4: (Color online) Experimental Qs value spectum 
from the selection of Fig[3] for the kinematic locus of 
?H(°Li, pt)*He. The relevant peak is the one at about 2.5 
MeV. 


examined. This gives a major constraint for the presence 
of the quasifree mechanism and the possible application 
of the THM. 


In order to extract the experimental momentum distri- 
bution of the spectator in the system where the Trojan 
horse particle b is at rest, narrow energy and angular 
windows should be selected. Since (do /dQ)TE is nearly 
constant in a narrow energy and 6... window, one can 
obtain the shape of the momentum distribution |W]? of 
the undetected particle directly from the three-body re- 
action yield divided by the kinematical factor KF, ac- 
cording to Eq.(5]. 


The obtained momentum distribution is reported in 
Figure [5] where it is compared with the theoretical pre- 
diction of the spectator momentum distribution, ob- 
tained using the Woods-Saxon potential with the stan- 
dard geometrical parameters. 


An evident distortion of the momentum distribution 
shows up and its measured full width at half maximum 
(FWHM) turns out to be around 23 MeV/c which is 
much smaller than the expected prediction of 72 MeV/c. 
This evidence was already observed in Ref. [13], where 
the width of the momentum distribution for the spectator 
inside the Trojan horse nucleus was studied as a function 
of the transferred momentum q from the projectile a to 
the center of mass of the final system B = C +c. In the 
present case, the value of q is about 133 MeV/c, and the 
width of the momentum distribution is about 23 MeV/c. 
It is in agreement with the trend of the curve that rep- 
resents the best fit to the function reported in Ref. 
Wrwam (4) = foll — exp(—4:/40)]. 

For the further analysis, the condition of |ps| < 
20MeV /c will be added to the above cuts to select the 
quasi-free events of the three-body reaction. 
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FIG. 5: (Color online) Experimental spectrum of momentum 
distribution for intercluster motion of a inside °Li (the blue 
points and the dotted green line for the fitting curve) com- 
pared with the theoretical calculation (the red line) 


D. S(E) factor and U. 


The last step is to extract the energy trend of the S(E) 
factor by means of the standard procedure of the THM 
after selecting the qusi-free three-body reaction events. 

Therefore, Eq.(5) and Eq.(6) are applied. The rele- 
vant two-body reaction cross section can be extracted 
from the measured three-body cross section with cut of 
selecting the quasi-free events from the three-body reac- 
tion. Then, the S(E) factor can be determined from the 
definition of S(E) = o(E)E exp(277), where the Sommer- 
feld parameter is n = Z1Z3€?/ (iv). In present work, only 
the s-wave (J=0) was considered for the energy range of 
Ecm = 0 — 400keV. 

The results for the d(d,p)t reaction in terms of the 
bare nucleus astrophysical Spare(E) factor are presented 
in Figure [6] (blue points) after normalization with direct 
data (red points) ) hid b2 The normalization was per- 
formed in the energy range of Ecm = 40 — 400keV, in 
which the electron screening effect is still neglectable. It 
should be pointed out that direct data suffer from the 
electron screening effect which does not affect the THM 
results. Thus, the S(E) extracted via THM was called 
Spare (E). A polynomial fit was then performed on the 
data giving Spare(0) = 56.7 + 2.0keV - b. 

The data from the present experiment (blue points) 
are compared with those from PRC-2013 of 9Li = 
(d + a) breakup in a previous THM experimental run 
(pink points) and those from APJ-2014 of ?He = 
(d + p) beakup experiment (green points). An overall 
agreement is present among both direct and indirect data 
sets, within the experimental errors. 

It should be pointed out that the errors in the present 
case are much smaller than in the case of PRC-2013 
using the same Trojan horse with higher beam energy 
and AE detectors for particle identification. 
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FIG. 6: (Color online) The S(E) factor obtained from THM 
measurement compared with direct data 


It is also in agreement, within the experimental er- 
rors, with the result using a different Trojan horse ?He 
1g. That is, data extracted via the THM applied to 9Li 
and ?He breakup are comparable among themselves. The 
Trojan horse particle invariance is confirmed in an addi- 
tional and independent case which was already observed 
in Ref. 3]. 

The lack of screening effects in the THM Spare(E) fac- 
tors gives the possibility to return the screening potential 
Ue from comparison with direct data using the following 
screening function with U, as free parameter. 


fia (E) = o5 (H2) /oy(E) ~ exp(rqUe/ E) (7) 


The result is shown in Figure[7] The red points are the 
direct data by Greife et al. (1995) [23]. The blue dashed 
line is the fitting curve of direct data (screened), and 
the green line is the fitting of THM data (unscreened) 
of the present work. Thus, we obtain a value of Ue = 
13.2 + 4.3eV, which is also in agreement with the one of 
Ref. [16] Ue = 13.4 + 0.6eV . 


V. SUMMARY 


A new investigation of the ?H(°Li, pt)*He reaction 
measurement to extrat information on the astrophysical 
Spare (E) factor and screening potential U, for the d(d, p)t 
reaction via the THM is presented in the present paper, 
shown in Table [I] 
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FIG. 7: (Color online) Fitting of S(E) to obtain Ue 


TABLE I: Comparison of d(d,p)t indirect study via THM. 


Work TH Eo EX So(E) Ue 
(MeV) (MeV) (keV - b) (eV) 
Present Li = (d--o) 9.5 0.089 56.74 2.0 13.2 + 4.3 
[5] ĉLi=(d+a) 14 0.866 75-21 - 
[6] *He=(d+p) 17 0.178 57.7418 13.4 € 0.6 


not influence the THM results. Thus, this work gives an 
updated test for the Trojan horse nucleus invariance. 
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